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Superconductivity occurs in systems that have a lot of negative charge: the highly negatively 
charged (CuC>2)~ planes in the cuprates, negatively charged (FeAs)~ planes in the iron arsenides, 
and negatively charged B~ planes in magnesium diboride. And, in the nearly filled (with negative 
electrons) bands of almost all superconductors, as evidenced by their positive Hall coefficient in the 
normal state. No explanation for this charge asymmetry is provided by the conventional theory of 
superconductivity, within which the sign of electric charge plays no role. Instead, the sign of the 
charge carriers plays a key role in the theory of hole superconductivity, according to which metals 
become superconducting because they are driven to expel negative charge (electrons) from their in- 
terior. This is why NIS tunneling spectra are asymmetric, with larger current for negatively biased 
samples. The theory also offers a compelling explanation of the Meissner effect: as electrons are 
expelled towards the surface in the presence of a magnetic field, the Lorentz force imparts them with 
azimuthal velocity, thus generating the surface Meissner current that screens the interior magnetic 
field. In type II superconductors, the Lorentz force acting on expelled electrons that don't reach 
the surface gives rise to the azimuthal velocity of the vortex currents. In the absence of applied 
magnetic field, expelled electrons still acquire azimuthal velocity, due to the spin-orbit interaction, 
in opposite direction for spin-up and spin-down electrons: the "Spin Meissner effect". This results 
in a macroscopic spin current flowing near the surface of superconductors in the absence of applied 
fields, of magnitude equal to the critical charge current (in appropriate units) . Charge expulsion also 
gives rise to an interior outward-pointing electric field and to excess negative charge near the sur- 
face. In strongly type II superconductors this physics should give rise to charge inhomogeneity and 
spin currents throughout the interior of the superconductor, to large sensitivity to (non-magnetic) 
disorder and to a strong tendency to phase separation. 

PACS numbers: 



I. INTRODUCTION 

In a perfectly conducting fluid, magnetic field lines 
move with the fluid (Alfven's theorem [l]]). Thus it is 
natural to infer that the expulsion of magnetic field 
lines from the interior of a metal making a transition 
to the superconducting state (Meissner effect |2|) is likely 
to be associated with radially outward motion of electric 
charge^. However, the conventional theory of super- 
conductivity predicts no radial motion of charge in the 
transition to superconductivity Rather, BCS-London 
theory postulates that a spontaneous azimuthal charge 
motion near the surface is generated (Meissner current) 
to compensate the magnetic field in the interior, without 
however explaining what is the driving force for such a 
motion nor how angular momentum is conserved 

Instead, the theory of hole superconductivity [(| pre- 
dicts that superconductors expel negative charge from 
their interior towards the surface as they enter the su- 
perconducting state @, HJ to lower their kinetic energy 
associated with quantum confinement [10j] , whether or 
not an external magnetic field is present. In the pres- 
ence of an external magnetic field, the radial charge mo- 
tion will 'drag' the magnetic field lines with it as in 
a classical plasma the magnetic Lorentz force on a 
radially moving charge acts in the azimuthal direction, 
and the deflected electron motion generates a magnetic 
field in direction opposite to the applied one. Thus, the 
theory offers a 'dynamical' explanation of the Meissner 



effect [llj. In the absence of applied magnetic field, the 
radial outflow of charge gives rise to a spontaneous spin 
current [l^, predicted to exist in the ground state of all 
superconductors |l3j . 

The hypothesis that the transition to superconduc- 
tivity is associated with expulsion of negative charge is 
supported by the observation that high temperature su- 
perconductivity appears to be favored in materials with 
substructures that have excess negative charge, namely 
the (CuC>2) = planes in the cuprates, the (FeAs)~ planes 
in the iron arsenides, and the B~ planes in magne- 
sium diboride. The observation that in high T c mate- 
rials normal-insulator-superconductor tunneling spectra 
are asymmetric, with larger current for a negatively bi- 
ased sample[ll [H, 

is further evidence that supercon- 
ductors have a tendency to expel negative charge. Fi- 
nally, superconducting materials almost alwaysj3xhibit a 
positive Hall coefficient in the normal state[l||, which 
indicates electronic bands almost full with negative elec- 
trons. 

Superconductivity arises in our theory when the Fermi 
level is close to the top of a band[I3|, and T c is enhanced 
when the ions are negatively charged[l8j]. Electron- 
hole asymmetric electronic polaron models describe the 
physics of pair formation [l9j, which in the low en- 
ergy sector reduce to a Hubbard model with correlated 
hopping [20]. The pair (bipolaron) is lighter than the sin- 
gle polaron in these models [2l[ because the hopping am- 
plitude increases with increasing local hole occupation 
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due to electron-hole asvmmet rvl22i , and this effect pro- 
motes pairing of hole carriers [2(| • The models are de- 
rived from basic atomic physics considerations of wide 
generality 22j , and the theory is proposed to apply to all 
superconducting materials [23[. 



II. THE TWO ROUTES TO THE MEISSNER 
EFFECT 

The fact that the Meissner effect is unexplained by 
the conventional theory is not generally recognized 11| . 
The canonical momentum of an electron with superfluid 
velocity v s is 



where r is in the plane perpendicular to the axis of the 
cylinder. Hence ?■ B = and fx (v x B) = — (r ■ v)B, 
and 

^{fxv) = — (f-v)B = - — —{^-r 2 )B (6) 

dV ' m e c ' 2m e c y dt ' y ' 

so that f x v = ~{e/2m e c)r 2 B, and the acquired az- 
imuthal velocity in moving out a distance r is 



V<f, 



2m„c 



-rB 



(7) 



Thus, to acquire the azimuthal speed Eq. (2) needed for 
the Meissner current requires the action of the Lorentz 
force over a radially outgoing motion to radius r = 2A^. 



P 



m e v s 



-A 



(1) 



with A the magnetic vector potential. In the BCS ground 
state the expectation value < p >~ 0, hence the super- 
fluid velocity is given by 



v s = A 

m e c 



e\i 



-B x n 



(2) 



The second equality in Eq. (2) applies to a cylindrical 
geometry, where h is the outward pointing normal of the 
lateral surface of the cylinder and B is the magnetic field 
along the axis of the cylinder. The London penetration 
depth Ai is given by[l] 
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47m, 



(3) 



where n s is the superfluid density. 

Eq. (2) embodies the Meissner effect [3]. However the 
BCS 'explanation' just outlined does not explain how the 
electrons are driven to acquire this velocity starting from 
a normal state where the average velocity is zero in the 
presence of a static magnetic field, nor how the mechan- 
ical angular momentum of the carriers of the Meissner 
current is compensated [llj . 



B. Meissner current from Faraday induction 

For an electron orbiting in a circular orbit of radius r, 
as an external magnetic field perpendicular to the orbit 
is applied, an azimuthal electric field E = (r/2c)dB/dt 
is generated by Faraday's law, and the velocity of the 
electron changes as 



dv 
~dt 



eE 
m, 



dB 



2m e c dt 



(8) 



so that for a magnetic field increasing from to B the 
extra velocity acquired is 



Av 



2m P c 



-B 



(9) 



which reduces to Eq. (2) if and only if the orbit has radius 
r = 2\l- Hence, the hallmark property of superconduc- 
tors, that the same Meissner current Eq. (2) results when 
a magnetic field is applied to an already superconducting 
metal or when a normal metal becomes superconducting 
in a pre-existent magnetic field, can be understood from 
the assumption that superconducting electrons reside in 
mesoscopic orbits of radius 2Al. A parallel reasoning 
leads to the development of a ground state spin current 
in the absence of applied fields [13[ (Spin Meissner effect) 
as we discuss in the following section. 



A. Meissner current from orbit expansion 

Consider an electron that moves radially outward from 
the axis of a cylinder in the presence of a magnetic field 
B parallel to the cylinder. The equation of motion is 



dv 
'~dl 



-vx B + F r 



(4) 



where the first term is the magnetic Lorentz force and 
the second term is a radial force arising from "quantum 
pressure" that drives the electron outward [10]. From Eq. 

(4), 



dv e ^ 

r x — = r x (v x B) 

dt m e c 



(5) 



III. THE TWO ROUTES TO THE SPIN 
MEISSNER EFFECT 

The superconducting condensate carries a charge den- 
sity en s . Thus, in a charge- neutral system the superfluid 
moves in a compensating background of positive charge 
density p — \e\n s . The interaction of the moving mag- 
netic moments of the electrons with the positive back- 
ground leads to a universal spin Meissner current with 
speed of magnitude [HI 

•«-5£r < 10 > 

as we will shows in what follows, which parallels the dis- 
cussion in the previous section. 
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A. Spin Meissner current from orbit expansion 



B. Spin Meissner current from Maxwell induction 



Consider a magnetic moment p along the z direction 
that moves radially outward with velocity v. It is equiv- 
alent to an electric dipole moment [24| 



P = - X fJL 

C 



(11) 



In the presence of the radial electric field of the cylinder 
E = 2ixpr = 2n\e\n s r (12) 
the electric dipole experiences a torque 

-* v -* v 

f = px E = (— X a) X E = — 2ir\e\n s f x (- x 0) (13) 
c c 

which causes a change in its angular momentum 



Hence 



dL d _ . ^ 

— — = m e — (r x v ) = t 
dt dV ' 



dv 2iren s _ ,v 
fx — = r x (- x p) 



dt m e 



(14) 



(15) 



Eq. (15) is identical to Eq. (5) if we define the 'effective' 
magnetic field 



B a = 2itn s p 



(16) 



and hence leads to the azimuthal velocity Eq. (7) with 
B a replacing B 



nen s 
vj, = rfi B 



(17) 



The same result for the azimuthal velocity Eq. (10) 
is obtained through a reasoning paralleling the second 
route to the Meissner current (Sect. IIB). 

Consider a magnetic moment p in an orbit of radius r, 
with p oriented perpendicular to the plane of the orbit. 
Assume a radial electric field grows from to a final value 
E. According to Ampere-Maxwell's law a magnetic field 
is induced by the varying electric field, satisfying 



V x B 



1 dE 
c~dt 



(21) 



which exerts an azimuthal force on the magnetic moment 



-* dv — ' -* 
F = m e —=V{P-B) 



(22) 



We have 



V(p- B) = {p-V)B + px (V x B) (23) 



In the geometry under consideration the first term in 
Eq. (23) is half the second term and points in opposite 
direction, so that 

F = m e f t =1mx(VxS) = 1| ( mx^ (24) 

and the azimuthal velocity acquired is 

1 



2m„c 



x E 



and for the electric field given by Eq. (12) 

Tr\e\n s _ _ 
V4, = n x r 



m e c 



(25) 



(26) 



with [lb = \e\h/2m e c the Bohr magneton, so that 
im s e 2 hr fir 



or 



V<j, 



2m 2 e c 2 8m e A| 



(18) 



where we have used Eq. (3) for the second equality in Eq. 
(18). The two electrons in a Cooper pair have opposite 
spin and orbit in opposite directions. The orbital angular 
momentum of each electron is 



8X 2 , 



(19) 



For r — 2\l, the azimuthal velocity Eq. (18) reduces to 
Eq. (10) and the orbital angular momentum is 



l= % -. 
2 



(20) 



For any other value of r, the orbital angular momentum 
Eq. (19) is not h/2. 



Trn s e 



hr 



8m e \ 2 L 



(27) 



in agreement with Eq. (18). Thus, just like for the 
Meissner effect, the same spin-current azimuthal speed 
Eq. (10) is obtained for a magnetic moment moving ra- 
dially outward a distance 2A l in the presence of a radial 
electric field Eq. (12) as for a magnetic moment orbiting 
at radius 2Al that is subject to a time-dependent radial 
electric field that grows from zero to its final value Eq. 

Note that the finding that the orbital angular momen- 
tum of the electron in the Cooper pair is H/2 (Eq. (20)) 
was not "built in" . Rather, it was derived (through two 
equivalent routes) from the hypothesis that the size of 
the orbit is 2Al, which in turn was inferred from the exis- 
tence of the Meissner effect, together with the reasonable 
assumption that the outgoing electron magnetic moment 
interacts with a positive background of equal charge den- 
sity as the charge of the superfluid (p = |e|n s ). 
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The magnitude of the magnetic field that will stop the 
spin current velocity of one of the spin orientations (the 
one that is parallel to the applied B) satisfies, according 
to Eqs. (2) and (10) 



eA; 



Am e X L 



-B x 



hence it is given by 



he 



m e c 



4^A| 



(28) 



(29) 



with <i>o = hc/2\e\ the flux quantum. Eq. (29) is essen- 
tially the lower critical field of a type II superconductor, 
H c i, that will drive the system normal^, and it coincides 
with B a , Eq. (16). The flux of the "stopping field" B s 
through the area of the orbit of radius 2Az, is precisely 
the flux quantum $o- 



IV. NEGATIVE CHARGE EXPULSION 

We have shown in the previous sections that expansion 
of the electronic orbits from a microscopic dimension to a 
mesoscopic radius 2Al describes the Meissner effect and 
predicts the Spin Meissner effect [131]. This expansion also 
gives rise to expulsion of negative charge from the interior 
of the superconductor towards the surface, as we discuss 
in what follows. 

In the normal state, electronic orbits carry zero orbital 
angular momentum on average, however each electron 
has an intrinsic (spin) angular momentum h/2. We can 
think of the spinning electron as a charge e orbiting at 
speed c in an orbit of radius given by the "quantum elec- 
tron radius" r q = S/(2m e c) . We have seen that as the 
orbit expands to radius 2A^ the orbital angular momen- 
tum acquired is also h/2 (Eq. 20). Thus the electron 
orbiting at radius 2Al with orbital angular momentum 
h/2 can be regarded as a magnified image of the spin- 
ning electron, with 'magnification factor' 2\i J /r q . Hence 
it is natural to conclude that the charge e will be corre- 
spondingly reduced by the factor r q /(2\L), so that the 
expelled negative charge density is 



P- 



'2X, 



(30) 



Eq. (30) implies that the spin current can be equivalently 
regarded as being carried by charge densities en s /2 or- 
biting at speed ±i>° or by charge densities p_ /2 orbiting 
at speed ±c (the same charge density orbits in each di- 
rection in both cases in the absence of applied magnetic 
field) . A similar result holds for the charge current as we 
discuss below. 

Indeed it can be shown 25] that the requirement that 
the theory be relativistically covariant leads to the con- 
clusion that a negative charge density p_ of magnitude 
given by Eq. (30) exists within a London penetration 



depth of the surface of superconductors. This negative 
charge was expelled from the interior of the superconduc- 
tor in the transition to superconductivity [HI, resulting in 
an interior positive charge density 



Pa 



2\ L 



P- 



(31) 



for a cylinder of radius R. The electric field generated 
by this internal positive charge density increases linearly 
with r, the distance to the cylinder axis, and reaches a 
maximum value 



E m = 2-npaR 



he 



(32) 



within a London penetration depth of the surface. Note 
that Eq. (32) is the same (in cgs units) as the "stopping" 
magnetic field Eq. (29) as well as the effective spin-orbit 
field Eq. (16). 

We can also understand the result Eq. (32) from the 
following heuristic argument. The expelled charge den- 
sity p- is related to E m by 



E,, 



P- 



AttXj 



(33) 



due to charge neutrality. The Meissner current in an 
applied magnetic field B has magnitude 



J 



n» em. 



Eq. (34) can be interpreted as the current created by the 
excess negative charge p_ moving at speed 



B 

E m 



(35) 



and suggests that superconductivity will be destroyed 
when v p _ reaches the speed of light. This will occur 
for B = E m , thus the value of the magnetic field that 
stops the spin current and destroys superconductivity, 
Eq. (29), yields the value of the electric field near the 
surface E m Eq. (32). 

From Eqs. (3), (10) and (32) it follows that the electro- 
static energy density due to the electric field E m equals 
the kinetic energy density of the spin current 



1 p2 

1 f 0\2 m 

—m e (v„)n s = — — . 

2 V aJ 8tt 



(36) 



The same relation exists, as is well known [ij, between the 
kinetic energy of the Meissner current and the magnetic 
energy density 



-m e (v s ) 2 n s 



EP_ 
8tt' 



(37) 



as can be seen from Eqs. (2) and (3). 
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V. ELECTRODYNAMICS OF CHARGE AND 
SPIN 

The foregoing considerations lead to the following four- 
dimensional equation in the charge sector [27] 



J - Jo = - 



4ttA| 



(A-Ao) 



with the current four-vector given by 
J = (J(f,t),icp{f,t)) 



(38) 



(39) 



with J the charge current and p the charge density, and 
the vector-potential four-vector given by 



A = (A(r,t),i<j>(r,t)) 



(40) 



with A the magnetic vector potential and <p the elec- 
tric potential, related by the Lorenz gauge condition 
DivA = 0, with Div = (V,d/d(ict)). The quantities 
with subindex are 



Jo = (0, icpo) 



A Q = (0,# Q (f)) 



(41a) 



(41b) 



with V 2 0o = — 47rpo and po determined by Eqs. (30)- 
(32). The spatial part of Eq. (38) is the ordinary Lon- 
don equation. From the fourth component of Eq. (38) 
and Maxwell's equations it follows that there exists an 
electrostatic field in the interior of superconductors that 
satisfies 



V 2 (E-E ) = -L(E-Eo) 



(42) 



with Eq the electrostatic field generated by the uniform 
charge density po- 

The charge current four-vector Eq. (39) is composed 
of the sum of spin current four-vectors 



J = J T + J L 
and the spin current four- vectors satisfy [2£ 



J a ~ Ja-o 



(Ac - A a0 ) 



with 



J a = (Jcicpc) 



Act = (A a ,i(f>c) 



(43) 



(44) 



(45a) 



(45b) 



J a = e(n s /2)vc is the component of the current of spin 
a and p a is the charge density with spin a. The spin 
potentials are given by|25j 



X L a x E(f,t) + A(r,t) 



(46a) 



Finally, the quantities with subindex are 
Jao = (Jco(r),icpco) 



Jao(r) 



cpo 



Po 



and 



Aco = (Acoirj^^aoir)) 



Aco(r) = X L a x E (r) 



0ao(r) = (t>o{r) 



(46b) 

(47a) 
(47b) 

(47c) 

(48a) 
(48b) 
(48c) 



These equations predict the existence of a spontaneous 
spin current flowing within a London penetration depth 
of the surface of the superconductor, with carrier den- 
sities (n s /2) and opposite spin flowing in each direction 
with speed Eq. (10), and a spontaneous electric field 
throughout the interior of the superconductor, of maxi- 
mum value given by Eq. (32) [251 ] . 



VI. ENERGETIC CONSIDERATIONS 

The kinetic energy of a pair of electrons of opposite 
spin due to the spin current is e p — m e (i>°) 2 . When the 
applied magnetic field approaches B s (Eq. (29)), one of 
the members of the pair doubles its speed and the other 
one comes to a stop, hence the kinetic ener gy of the pair 
doubles, at which point the pair breaks up[lj]. We con- 
clude from this argument that the condensation energy 
of the pair is m e (?4J) 2 and hence that the condensation 
energy per electron is 



1 , 
—mAv 

2 y ' 



°) 2 



(49) 



which equals the electrostatic energy cost per electron 
due to the internal electric field, Eq. (36). This implies 
that each electron lowers its energy in entering the con- 
densate by 



2e r = 



4m K 



(50) 



with qo = 1/2 Al- The system gives back half of this gain 
right away in the electrostatic energy cost Eq. (36), and 
the other half when the applied magnetic field destroys 
superconductivity. 

What is the physical origin of this energy lowering? 
Note that v can be written as 



1 



7PBB a 



2m e c 2 



\S.(%xE)\ 



(51) 
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with S = h/2 and E the radial electric field Eq. (12) 
at radius r = 2Al, normal to v^. The second form of 
Eq. (51) is the usual spin-orbit energy including the cor- 
rection for Thomas precession (28j. From this we con- 
clude that the condensation energy of the superconductor 
originates in the spin-orbit energy lowering arising from 
the interaction of the spin current of the condensate ( of 
charge density en s and velocity Eq. (10)) with the com- 
pensating positive background charge density \e\n s . 

It is also interesting to note that v determines the frac- 
tion of the superfluid charge density (en s ) that is expelled 
(as suggested in [7[), through the relations 



P- 



en s (- 



e\E m \ L ' 



en s (- 



m e c^ 



1/2 



(52) 



where eE rn \h is the electrostatic energy difference be- 
tween an electron at the center and at radius 2Al of a 
cylinder with charge density (The second expres- 

sion in Eq. (52) was also found in ref.[26( through dif- 
ferent arguments), v is also related to this electrostatic 
energy difference through 



{eE m \ L ) 2 



(53) 



The parameter v represents a change in the chemical 
potential between normal and superconducting states @, 
8], which according to the theory of hole superconduc- 
tivity is related to the slope of the energy-dependent gap 
function A k bydHUHfl 



1 d 

2dT k 



(A fc ) 2 



(54) 



From Eqs. (50), (54) and (10) and using a free-electron 
dispersion relation e k = h 2 k 2 /2m e we find that the 
energy-dependent gap function is given by 



A k = 



h 2 q n k 
2m e 



= V2 



(55) 





(a) 



(b) 



FIG. 1: Schematic depiction of spin-split bands (a) and as- 
sociated real-space picture (b). (a) The spin up electrons 
(depicted as circles) occupy predominantly negative k-s and 
the spin down electrons (depicted as crosses) positive k's. The 
difference in energy between the highest occupied k-state for 
positive and negative k for a given spin orientation is twice 
the energy gap. (b) The excess negative charge near the sur- 
face of the superconductor flows predominantly clockwise or 
counterclockwise depending on its spin orientation. 



with the minimum quasiparticle gap Ao = y/2fMU — v 2 . 
For the second equality in Eq. (57) we used Eq. (55) . Eq. 
(57) shows that indeed v is the change in the chemical po- 
tential in going from the normal to the superconducting 
state, as anticipated in the "correlated hopping" model 
of hole superconductivity!?], 

Furthermore note that a spin current with speed 
Eq. (10) can be represented by the energy- wavevector 
relation 321 



C-ko 



-k,- 



2m e K 2 ' 



(58) 



with a = ±1, where qo is a vector in the direction of 
the spin current flow of magnitude go = 1/2Al (since 
h-^de^/dk - de kl /dk)/2 = hq /(2m e ) = «°). Eqs. 
(58) and (55) then imply that 



Both the parameter v and the gap at the Fermi energy 
A kp can be expressed in terms of the slope of the gap 
function at the Fermi energy, to, that determines the 
tunneling asymmetry 



5Aa 

de k 



)e 



go 

2k F 



2meF 



2m tp = mA k 



(56a) 



(56b) 



(56c) 



which illustrates that a sloped gap function is a necessary 
condition for superconductivity [Til. [30L l3l| . 

The quasiparticle energy is given within BCS theory 

by 



E 2 = ( Cfc - //) + A 



(e fc - /i + u) 2 + A' 2 (57) 



A fc 



(59) 



for k || <fo- In other words, the superconducting energy 
gap is due to 'spin splitting' [33|. Note also that the 
'depairing' speed which will cause the spin current to 
stop and the pairs to break is given by Eq. (10), which 
can be written as 



(60) 



hkp 



in terms of the gap Eq. (55) at k = k F - Eq. (60) 
for the critical speed is identical to what is obtained in 
conventional BCS theory Q, even though the energy gap 
expression Eq. (55) was obtained through an entirely 
independent argument. Fig. 1 shows schematically the 
spin-split bands in momentum space and the associated 
charge configuration and spin current in real space. 



(a) Type I superconductor 



(b) Type II superconductor 



FIG. 2: Schematic depiction of a Cooper pair in a type I (a) 
and type II (b) superconductor. The vertical arrows denote 
the direction of the electron magnetic moment, and the hor- 
izontal arrows the orbiting direction. The radius of the orbit 
of each electron is 2Al, and the distance between the centers 
of the orbits is £. In type II materials ((b)) with £ < 2Al 
a normal vortex core of diameter £ can be enclosed by both 
orbits of the same Cooper pair. 



Finally, note that the condensation energy per electron 
Eq. (49) can be written, using Eq. (55) and the free- 
electron dispersion relation, as 



4e_F 



(61) 



Eq. (61) is consistent with the BCS expression for the 
condensation energy per unit volume Q 



SU = -N{0)A 2 



(62) 



with iV(0) the density of states per spin at the Fermi 
energy, if we take -/V(O) = n s /2eF appropriate to a two- 
dimensional free-electron system. 

Van der Marel(34j and Khomskii[35| have pointed out 
that quite generally in conventional BCS theory a shift 
in the chemical potential is predicted upon entering the 
superconducting state, of magnitude A 2 /Aep- In our case 
the shift in the chemical potential, v, is twice as large 
(Eqs. (50) and (61)), and it is directly related to the 
slope of the gap function and to the existence of a spin 
current and of negative charge expulsion. 



VII. TYPE I VERSUS TYPE II MATERIALS 
AND CHARGE INHOMOGENEITY 

The above energetic considerations apply close to the 
crossover between type I and type II behavior, where 
H c ~ H cl and Xl ~ £, where £ is the coherence length 
which is also the average distance between members of 
a Cooper pairQ. In extreme type I materials, pairs are 
broken well before the charge speed reaches Eq. (10), 
<?o ~ 1/? << 1/2 Xl in Eq. (55) and E m = H c rather 
than Eq. (32). We can understand the crossover between 
type I and type II behavior geometrically, as shown in 
Fig. 2. £ is the distance between the centers of the 2Al 
orbits of up and down spin electrons in a Cooper pair. A 
vortex core has diameter £, and this normal region can 
be enclosed by the orbits of both members of the same 
Cooper pair if and only if £ < 2A^. 




FIG. 3: Schematic depiction of a superconductor with strong 
disorder in the absence of applied magnetic field. Defects, 
grain boundaries, vacancies, etc. will result in patches of 
normal regions (hatched areas) surrounded by spin currents 
(dashed lines, with arrows pointing in the direction of flow 
of electrons with magnetic moment pointing out of the pa- 
per) and excess negative charge density (gray areas). The 
figure also shows the excess negative charge and spin current 
near the surface. If the system is cooled in the presence of a 
magnetic field, magnetic flux will be trapped in the hatched 
regions and a charge current will flow around those regions to- 
gether with the depicted spin currents. The smallest normal 
regions have diameter of a coherence length. 



Why does the vortex core have to be enclosed by the 
orbits of both members of the same Cooper pair? The 
phase change for each electron in going around a loop 
enclosing a vortex is 



li\H = 4> m e v ■ dl + -(f>B 



(63) 



where c/)b is the enclosed magnetic flux. For each mem- 
ber of the Cooper pair this phase change is 7r, corre- 
sponding to its angular momentum h/2, and thus 4>b — 
(h/2)c/e — 4>q (assuming the integration loop is through 
a path where v = 0), thus providing a new rationale for 
the factor of 2 in the flux quantum (f>Q. If only one mem- 
ber of the Cooper pair were to enclose the vortex the 
wave function for the pair would not be single- valued. 

The superconductor expels negative charge towards 
the surface and towards any interior normal regions, thus 
in the flux phase there will be excess negative charge in 
and around the vortex cores. Furthermore our model 
has a large sensitivity to disorder arising from the slope 
of the gap function, as discussed in [36[. In the presence 
of local potential variations due to impurities, vacancies, 
etc, the gap can be sharply reduced or vanish altogether 
giving rise to normal regions. There will be excess nega- 
tive charge in and around those normal regions expelled 
from the superconducting regions, and a spin current will 
circulate around the interior normal regions as shown 
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schematically in Fig. 2. An applied magnetic field will 
concentrate in these normal or weakly superconducting 
regions and the spin current around them will acquire a 
charge current component. 

Finally, the model has a strong tendency to phase sep- 
aration in the regime where the carrier (hole) concentra- 
tion is small (37}. This arises from the bandwidth depen- 
dence on the hole concentration[38j: the bandwidth in- 
creases with increasing hole concentration, favoring seg- 
regation into hole-rich and hole-poor regions. This ten- 
dency is also enhanced in the presence of disorder. 

VIII. SUMMARY AND DISCUSSION 

The theory discussed here proposes that superconduc- 
tivity arises from the fundamental charge asymmetry of 
matter [3^|. It was motivated by the discovery of high 
T c cuprates 20] but it was clear from the outset that if 
valid it would apply to all superconducting materials [iol]. 
As the theory progressed it became increasingly appar- 
ent that it led to a radical departure from the conven- 



tional BCS theory [1, [2l|. Later it became clear that 
even London electrodynamics had to be modified [13] ■ Fi- 
nallyit became clear that spin-orbit coupling plays a key 
rolejl3l Hfl . The fact that at the end of this road one is 
led, unexpectedly, to a far more compelling (in our view) 
explanation of the most fundamental property of su- 
perconductors, the Meissner effect, than the conventional 
theory proposes, is in our view a strong argument for its 
validity. Ultimately of course confirmation or refutation 
of the theory discussed here will come from experiment. 

As the theory advanced we have gained increased un- 
derstanding of what we knew we didn't know, however 
we have also become aware of issues that previously had 
been 'unknown unknowns' [4lj. Perhaps the most impor- 
tant of these issues, that remains to be understood, is the 
proper inclusion of the key role of the Dirac sea[42j|. 
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